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G  shear  modulus  of  the  matrix 


ABSTRACT 


A  number  of  authors  have  employed  shear  lag  theory 
to  find  the  stress  distribution  near  cracks  in  a  uni- 
ax  tally  reinforced  composite.  Practical  application  of 
these  equations  has  been  hampered  by  the  fact  that  they 
include  a  parameter  h/d  the  magnitude  of  which  is  un¬ 
known.  This  parameter  is  evaluated  for  composites  in 
which  the  fibers  form  a  rectangular  array. 

NOMENCLATURE 

d  distance  between  rivet  lines  or  distance  between 
rod  centers 


d^  distance 

d^  distance 

h  kvlate  or 

j  jr  of  the 

dc.h  d 

h  h  , 

d"  d  of  th€ 

c.v 


between  next-to-nearest  fiber  pair 
between  nearest  fiber  pair 
shell  thickness 

composite  In  the  horizontal  direction 

composite  in  the  vertical  direction 


j  current  density  in  electrolyte 

Jn  normal  component  of  the  current  density 

r  radius  of  a  fiber  or  a  rod 

o 

r^  distance  from  a  negative  line  force 

r^  distance  from  a  positive  line  force 

w  axial  displacement  of  the  n'th  fiber 
n 

z  axial  coordinate 

Ft  force  in  the  I'th  stiffener  or  fiber  in  a  composite 


Gh  effective  shear  stiffness  of  the  composite 
I  current  generated  by  the  power  supply 
1^  current  in  the  i'th  rod 


Pe  resistivity  of  the  electrolyte 
T  shear  stress  In  matrix 

♦  electrical  potential  in  rod  or  electrolyte 
INTRODUCTION 

The  analysis  of  failure  of  initially  crack-free 
continuous- fiber  composites  is  complicated  because 
such  failures  are  not  governed  by  weakest  link  cons id 
eratlons.  Instead  they  are  a  result  of  damage  accumu 
la t ion.  To  account  for  crack  growth  we  must  know  the 
microstresses  in  the  neighborhood  of  the  crack  tip. 

At  present  such  understanding  Is  limited  to  unidirec- 
tlonally  reinforced  composites  and.  even  there,  is 
somewhat  fragmentary. 


Most  attempts  to  find  stress  distributions  near 
the  crack  tip  are  based  on  shear  lag  theory.  This 
approach  simplifies  elastic  theory  by  assuming  that 
all  displacements  are  parallel  to  the  fibers  and  that 
the  matrix  between  the  fibers  carries  only  shear 
stress.  The  fundamental  equation  of  shear  lag  theory 
originally  developed  for  the  stiffened  plates  and 
shells  employed  In  aircraft,  takes  the  form 


Hi 

dz 


-  Gh  — 


*i-l 


+  2w^  - 


i+1 


(IV 


By  analogy,  the  corresponding  equations  for 
composites  were  assumed  by  Hedgepeth  (1,2)  and  later 
by  many  others  to  take  the  same  form  in  the  case  of  a 
linear  array,  and  in  the  case  of  a  square  array  of 
fibers  to  take  the  form 


ds 
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♦Adjunct  Professor,  Fallow  A1AA 
♦♦Research  Assistant,  Student  Member  AlAA 


This  would  be  correct  for  a  square  array  if  it  had 
a  very  low  fiber  volume  ratio  and  the  matrix  was  re¬ 
placed  by  thin  plates  as  shown  in  Figure  1.  But  for 
high  fiber  volume  ratio  and/or  conventional  matrices 
it  is  quite  unclear  how  to  estimate  h/d  or  whether 
the  assumption  that  only  nearest  neighbors  interact 
i-  .%  iood  one. 

iotiu  work  has  already  been  done  on  these 
questions.  Batdorf  (3)  gave  an  exact  solution  for 
the  interaction  between  two  infinitely  long  rigid 
rods  immersed  in  an  Infinite  elastic  matrix  with 
displacements  t  The  result  in  this  case  was  that 


C  ~  •  2w 
d  o 


-/MU  i  vl-(2r  /d)l/(2r  /<*)  ] 
o  o 


Here  rG  is  the  rod  radius  and  d  is  the  distance  between 
rod  centers.  It  follows  that  h/d  is  a  function  of  r0/d 
varying  from  0  for  r0  -  0  to  «  for  rQ  -  d/2.  Unfortu¬ 
nately  the  technique  employed  cannot  be  extended  to 
more  than  two  fibers.  This  is  because  whereas  the  dis¬ 
placement  contours  in  the  matrix  surrounding  two  rods 
of  circular  cross-section  are  circular  cylinders,  for 
more  than  two  rods  they  are  not. 

The  interaction  between  nearest  neighbors  for  the 
case  of  a  square  ^rray  was  found  by  Batdorf  and 
Ghaffarian  (4)  using  an  electric  analogue  (5) .  The 
present  papers  extends  this  approach  to  the  determina¬ 
tion  of  interaction  between  nearest  neighbors  in  a 
rectangular  array.  To  supplement  the  solution  found 
experimentally  employing  the  electric  analogue,  an 
analytical  approach  valid  for  very  snail  fiber  volume 
ratio  is  also  presented.  The  two  solutions  are  found 
to  be  in  good  agreement  in  the  limiting  case  where  the 
analytical  solution  is  valid. 


Analytical  Solution  for  Rectangular  Array 

For  the  case  of  a  rectangular  arrav,  two  parameters 
have  to  be  taken  into  account  at  the  same  time;  namelv, 
for  the  two  nearest  neighbors  and  4  for  the  next-to- 
nearest  neighbors.  The  analysis  for  such  an  array  can  he 
carried  out  by  a  consideration  of  the  load  transfer  du*  t 
a  set  of  equal  but  opposite  forces  directed  *.  te’  r. 
neighbors  as  in  Figure  2  and  r<>  twe  next  -r  -''em--* 
neighbors  as  in  Figure  3.  The  jnalysis  is  valid  11*  •: 
limiting  case  of  small  rn/d. 
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FIG.  1  HYPOTHETICAL  MATERIAL  OBEYING  CONVENTIOV.L 
cqr&n  Mr  EQUATIONS 


FIG.  3  LOADS  APPLIED  **'0  NEXT-TO-M.APd 
NETGHBOk  PAIR  -  CASF  B 

In  Case  A,  the  shear  lag  equation  becomes 

Fa  '  Gl(?)c,v  (2  G  ‘  Wf  ’  V  1 

(?)  .  (?  «  -  v.  -  w  i  1 
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-  G  (7)  <«,  *  V 

d  i,v  8  C 

where  subscript  c  denotes  the  value  for  the  2-fiber 
interaction  in  a  rectangular  array  composite,  and  i 
denotes  the  value  for  a  pair  of  isolated  fibers.  For 
sufficiently  small  r0/d2, 

w.  ’  -  wc  *  fe  ln  (Vro  "l)  (6 


,K  n  ^Viw2  x  , 
(d>  (1  ’  MdTTr -T)}  + 

c,v  1  o 


/  h  \  /  i  c  /2 _  . 

d0  (1,5  '  *n(d  /r  -1)  } 

c,  n  1  O 


The  displacements  of  the  other  fibers  can  be  obtained  via 

wi  "  fnG  in  ^rl/r2^  ■  1  **  a  or  c 

■where  the  field  point  is  the  center  of  the  fiber  whose 
displacement  is  being  calculated.  Thus 


b  2  WG 


Vl+l/a*  — 


in  Yl+l/a2 


in[(l  ♦  Vl-(2r  /d.  )  2)  /  (  2r  /d.  )  ] 
o  l  0  1 


Solving  (13)  and  (14)  simultaneously,  we  obtain 

,h.  _  A  -  (B/F)  C 

V  .  D  -  C(E/F) 
c  ,h 


8  -  Cj)  (I) 
- £.s.'ti _ 


(by  symmetry) 


1  -(2r  /d,)  )/(2r  /d,)J 


e  a  in(d-/r  -1) 
Z  o 


Thus,  (5a)  can  be  written  as 


F;-G'<7>c,v  <m/2  -  ln(lh  -1)>  + 

t  o 


,tu  /,  in V 14-1 /a2  x 

(J}  .  (1  ~  S"(d,7T~i))1  (2V 

C  ,n  Z  o 


=■  GC^)  (u  -  w  ) 

1  i,v  8  C 

Employing  (4)  for  £  for  a  pair  of  Isolated  fibers, 
(*>  can  be  replaced  by 
di,v 


in((l±VL-(2r  /d0)*)/(2r  /dj  ] 
o  Z  o  i 


Again  by  noting  that  wa  ■  -we,  the  above  equality  be¬ 
comes 


(t)  (1.3  - 


/h,  n  tn  Vl+l/tt*  ) 

'd2  .  u  "  tn(d,/r  -1) 

c ,  n  Z  o 


in((li>^-(2r  /d.)2)/(2r  /d.)]  (13) 

o  i  o  i 

The  above  equation  is  a  linear  algebraic  equation  In 
two  unknowns  (h/d)c  v  and  (h/d)c  ^ .  A  second  equation 
can  be  found  by  applying  a  similar  approach  to  Case  B. 


!n[(l  tV  1  -( 2«r  /  d. ) 2)  /  (2ar  Id.)) 

o  i  o  I 


C  -  1  - 


in  /1+a7 

in(d. /r-1) 
1  o 


D  -  1.5 


tn(dj/ro-Tj 


p  -  i  *n»A~+l/aT 

“  *  *  “  «  #?  nr  « \ 


in(d. /ar  -1 ) 
1  o 


F  *  1.5  - 


in(d./ar  -1) 
1  o 


Theory  Underlying  Electric  Analogue 

If  a  set  of  n  infinitely  long  rigid  rods 
Immersed  in  an  infinite  elastic  matrix  are 
displaced  axially  from  their  equilibrium 
positions  by  distances  w^,  w  . . .w  ,  the 
distorted  shape  of  every  cross-secP ion  will 
be  the  same.  Therefore  there  will  be  no 
direct  stress  In  the  matrix,  and  the  only 
shear  stresses  will  be  txz  *  tzx  and  XyZ  * 

TyZ .  At  the  Interface  between  the  iftn  rod 
and  the  matrix  the  equilibrium  equation  is 


dFt  L 

d T  ’  f  (T«  dy  +  Tzy  dx) 


where  0  signifies  Integration  around  the  circum 
ferenci  of  the  i’th  rod.  In  the  matrix,  equi¬ 
librium  In  the  *-direction  requires  that 


Dx  3x 

_i*  +  _ £2 

3x  3y 


(18) 


Since  the  only  displacement  is  w(x,y),  application  of 
Hooke's  lav  leads  to 


3xa 


+  32w 


3yl 

In  the  matrix,  and 
dF 


ati  l 

—  »  G  <p  >3v  ,  3v  . 

d*  Tt  <37  d*  +  dy) 


(19) 


(20) 


at  the  interface  with  the  i'th  fiber. 


Next  consider  a  geometrically  similar  array  of  n 
i  erfectly  conducting  rods  Immersed  In  an  electrolyte  of 
resistivity  pe.  Let  the  potentials  applied  to  the 
rods  be  $2  •  *  •  $n«  The  current  flowing  from  the 

electrolyte  into  a  unit  length  of  the  i'th  rod  is  then 

j>  0y  <*x  +  it  dy) 
i 

In  the  electrolyte  Ohm's  law 


-  lit 

dz 

leads  to 


(21) 


°e  ^  *  grad  # 

and  for  steady  currents 

7  •  j  -  0 

Thus  in  the  electrolyte 


+  i!i 

3x2  3y2 
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and  at  the  periphery  of  a  rod 


(22) 

(21) 


(24) 


dz 


(25) 


Comparing  (19)  and  (20)  with  (24)  and  (25)  we  see  that  $ 
is  analogous  to  w,  I  is  analogous  to  F,  and  pe  is 
analogous  to  1/G. 


The  analogy  provides  a  convenient  means  for  an 
experimental  determination  of  the  Interaction  between 
fibers.  If  we  assume  the  validity  of  (2),  we  can  write 


7{-Vj. 


(26) 


If  vc  measure  the  total  current  per  unit  length 
in  the  i'th  rod  when  known  potentials  are  applied  to 
that  rod  and  its  four  nearest  neighbors,  we  can  use 
(26)  to  determine  h/pcd.  By  measuring  pc  we  can 
find  h/d  for  a  square  array.  This  was  done  in  (4). 


The  case  of  a  rectangular  array  is  a  hit  more 
complicated.  The  square.  array  leads  to  (2b),  a 
linear  algebraic  equaticn  in  one  unknown,  h/d.  As 
noted  in  the  preceding  section,  in  a  rectangular 
array  there  are  two  unkrowns,  h/d  for  the  two 
nearest  neighbors  and  h/d  for  the  next-to-nearest 
neighbors.  Thus  two  experiments  must  be  carried 
out,  leading  to  two  simultaneous  equations  in  two 
unknowns . 


Experimental  Setup 

The  structure  of  tie  composite  was  simulated 
in  the  electric  analogue  by  conducting  rods  re- 
o lacing  the  fibers  and  a  weak  electrolyte  replacing 
the  matrix.  A  transparent  plastic  tank  was  used  to 
contain  the  fluid  and  the  rods. 


The  experiment  was  carried  out  first  by  in¬ 
serting  one  end  of  the  rods  into  a  plexiglass  panel 
with  holes  that  match  the  size  of  the  rods  and  with 
spacing  that  meets  the  rG/d  and  ot  under  investigation. 
This  set  up  was  then  loaded  into  the  electrolyte. 
Alignment  of  the  electrodes  was  provided  by  inserting 
the  other  ends  into  a  similar  panel  that  can  be 
mounted  at  any  convenient  level  above  the  elect rolvte. 

Two  experiments  were  performed  for  each  i  and 
each  r0/d.  The  first  one  involved  the  application  of 
a  constant  AC  potential  across  rods  a  and  b  as  in 
Figure  4.  The  two  outermost  rods  were  connected 
and  thus  d  was  at  the  same  potential  as  b,  while  c 
and  e  were  left  free  to  adopt  the  local  potential 
of  the  electrolyte. 


The  governing  equation  can  be  written  as 


« . 

dz 


i'fy  <V*ae)+(?>  <*.b+*.d)1 

e  c,v  c,h 


(?7) 


I  was  obtained  by  measuring  the  current  generated  bv 
the  power  supply,  z  was  the  submerged  length  of  the 
rods.  and  were  obtained  by  measuring  the 

potential  difference  between  a  and  c,  and  a  and  b , 
respectively.  By  symmetry,  t»ae  -  \xc. 


FIG.  U  POTENTIALS  APPLIED  TO  NEXT-TO-NFARF.ST 
NEIGHBOR  PAIR 
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In  order  to  obtain  a  second  equation  relating  the 
two  unknowns,  (Jj)  and  (il)  ,  a  constant  potential  vas 
d  £  y  d  C  Jl 

next  applied  acrosi  rods  a  ana  e  with  e  connected  to  c 
as  In  Figure  5.  The  preceding  governing  equation  again 
applies-  Similarly,  It  and  4>ab  were  measured  for 
this  case.  Hy  solving  the  two  linear  algebraic 
equations  simultaneously,  values  for  (*D  and 
<t)  were  obtained.  d  c,v 
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FIG.  11  h/d  FOR  RECTANGULAR  ARRAY 


Figure  6  gives  h/d  as  a  function  of  r0/d}  for  a  square 
array ,  i.e.,  a  *  1.  This  case  was  treated  earlier  in 
Reference  4.  The  other  figures  treat  the  cases  ~t  -  l.J 
to  a  »  1.5.  For  a  >  1  the  array  Is  rectangular  and 
therefore  there  are  two  curves,  one  representing  (h/d) 
for  the  vertical  interaction  and  the  other  for  the 
horizontal  interaction.  In  each  case  d  in  the  abscissa 
was  chosen  to  be  the  longer  distance  d x ;  this  was  of 
course  an  arbitrary  decision.  In  the  case  of  the 
ordinates  d  is  dy  for  the  horizontal  interaction  and 
d2  for  the  vertical  interaction. 

In  every  case  h/d  approaches  zero  as  r0/di 
approaches  zero,  and  in  this  region  the  theoretical 
solution  is  valid.  It  is  therefore  used  to  estimate 
the  value  of  h/d  at  such  small  values  of  v0/t*y  as  to 
make  an  experimental  determination  difficult. 

As  r 0/<*i  +  0.5,  the  horizontal  separation  be¬ 
tween  two  fibers  vanishes  and  (h/d)c  h  -*•  *».  This 
asymptote  can  be  used  to  extrapolate  ’the  experi¬ 
mental  data  to  higher  values  of  r0/d}  than  can 
readily  be  Investigated  experimentally.  As  r0/d2  * 

9.5  the  value  of  (h/d)c  v  *♦  00  because  then  the 
vertical  separation  shrinks  to  zero.  In  this  case 
rD/di  *  r0/ad2  •  0.5/a,  which  serves  similarly 
as  an  Asymptote. 

Discussion 

In  Hedgepeth's  papers  (1.2)  no  attempt  was  made 
to  evaluate  h/d  for  a  particular  composite,  and  no 
technique  was  proposed  for  doing  so.  These  paper? 
evaluated  the  stress  concentration  factors  in 
fibers  adjacent  to  a  crack,  and  it  turns  out  that 
the  stress  concentration  factor  does  not  depend  on 
the  value  of  h/d.  However  to  find  the  stresses  at 
which  crack  extension  occurs  it  is  necessary  to  find 
lengths  of  the  overloaded  fiber  elements,  and  these 
do  depend  on  h/d.  In  the  absence  of  a  theoretical 
technique  for  accomplishing  this,  various  assump¬ 
tions  have  been  made.  The  most  common  has  been  to 


I 


Assume  that  for  the  vortical  interaction  (interaction 
between  cloaeat  neighbors)  d  -  d2  ahd  h  -  dj ,  while 
for  the  horltontal  interaction  d  -  dj.  and  h  ■  d2» 

Refs.  (fcS) •  These  assumptions  lead  to  the  results 

(h/d)  -  a  (28a' 

c*v 

(h/d),.  .  •  l/o  (28b> 

c,h 

a  result  independent  of  r0/d.  Comparing  these  results 
with  Figures  6*11  we  see  that  they  are  somewhere 
near  right  for  rp/d^  -  0.2,  but  are  too  high  for 
smaller  values  of  rQ/dx  and  too  low  for  higher  values 
of  r0/dx*  We  note  that  when  r0/di  Is  0.2  the  fiber 
volume  ratio  is  about  a/8(  a  lower  value  than  what  is 
usually  encountered  in  practical  applications.  For  a 
square  array  having  a  fiber  volume  ratio  0,5, 
r0/d  3  0.4,  in  which  case  the  value  we  obtain  (see 
Figure  6)  is  about  2.5  times  that  given  by  (28). 

Concluding  Remarks 

During  the  last  three  decades  shear  lag  has  been 
widely  used  as  a  technique  for  finding  the  stresses 
in  damaged  unlaxially  reinforced  composites.  The 
applicable  equations  contain  a  parameter  representing 
the  shear  Interaction  between  neighboring  fibers. 
Somewhat  surprisingly,  up  to  now  little  has  been 
known  about  how  to  evaluate  this  term  for  a 
particular  composite. 

It  is  shown  that  for  small  values  of  the  fiber 
diameter  to  fiber  separation  this  parameter  can  be 
evaluated  theoretically.  For  larger  values  it  can 
be  evaluated  experimentally  employing  an  electric 
analogue. By  combining  the  two  techniques,  the  para¬ 
meter  is  evaluated  over  the  entire  range  of  diameter 
to  separation  ratio  for  composites  in  which  the  fibers 
form  a  rectangular  array.  The  results  contained  in 
Figures  6-11  cover  arrays  ranging  from  square  to 
rectangles  having  a  side  ratio  of  1.5. 
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